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Ferromagnetic coupling and magnetic anisotropy in molecular Ni(II) squares 
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We investigated the magnetic properties of two isostructural Ni(II) metal complexes [Ni 4 Lg] and 
[NULI]- In each molecule the four Ni(II) centers form almost perfect regular squares. Magnetic 
coupling and anisotropy of single crystals were examined by magnetization measurements and in 
particular by high-field torque magnetometry at low temperatures. The data were analyzed in terms 
of an effective spin Hamiltonian appropriate for Ni(II) centers. For both compounds, we found a weak 
intramolecular ferromagnetic coupling of the four Ni(II) spins and sizable single-ion anisotropies of 
the easy-axis type. The coupling strengths are roughly identical for both compounds, whereas the 
zero-field-splitting parameters are significantly different. Possible reasons for this observation are 
discussed. 

PACS numbers: 33.15.Kr, 71.70.-d, 75.10.Jm, 75.30.Et, 



I. INTRODUCTION 

Modern inorganic chemistry provides arrangements of 
magnetic metal ions in highly symmetric geometries, 
where the metal centers are separated by organic ligands. 
Compounds with a topologically simple arrangement of 
their metal centers like grids and rings have been inves- 
tigated extensively. 0, 0, 0i 01 IE El Due to their high 
symmetry it is possible to experimentally determine the 
magnetic parameters of large and rather complex sys- 
tems. As intermolecular effects are negligible in these 
systems, they actually form perfect models to explore 
finite size spin systems. In a number of cases, sev- 
eral related species of a compound class are available, al- 
lowing to observe correlations between crystallographic 
structure and magnetic properties. |fj, |8( Investigation 
of such correlations is of key help for understanding e.g. 
coupling mechanisms in detail. [!j 

Recently, the highly symmetric Ni(II) square [Ni 4 Lg] 
with L fc = C5H4N-CON-CN4-C2H5 attracted consider- 
able interest. 0, |n| Preliminary magnetic studies of 
powder samples revealed a sizeable ferromagnetic cou- 
pling of the four Ni(II) metal centers within a molecule. 
This system is thus one of the very rare examples of a 
ferromagnetic Ni(II) complex. A second species, [NljXg] 
with L c =^H4N-CSN-CN 4 -C 2 H5, could be also syn- 
thesized. [Ill 03 The two ligands used differ only at 
one position, i.e the oxygen in L h is replaced by a sul- 
fur in L c . Interestingly, the two Ni 4 squares are not only 
isostructural. Their magnetically relevant geometrical di- 
mensions actually differ by less than 3.5%. This leads to 
a new situation: Differences in magnetic properties stem 
predominately from different electronic properties of oxy- 
gen and sulfur and not from geometrical distinctions. It 
is the purpose of this work to investigate the magnetism, 
especially the magnetic coupling and anisotropy, of these 
two compounds in detail. 

When focusing on anisotropic properties, several ex- 



TABLE I: Comparison of selected distances and bonding an- 
gles of the coordination sphere for the compounds [Ni 4 Lg] and 
[Ni 4 L§]. Considering the nickel centers at the upper right cor- 
ners in Figs.Q^b) and|TJc). respectively, the N atoms of their 
coordination spheres were numbered clockwise from Nl to N5 
starting with oxygen/sulfur. 









[Ni 4 L|] 


[Ni 4 L|] 


Ni- 


O/S 




2.032(2) A 


2.051(3) A 


Ni - 


Nl 




2.081(3) A 


2.073(4) A 


Ni - 


N2 




2.076(3) A 


2.072(4) A 


Ni - 


N3 




2.111(3) A 


2.074(4) A 


Ni - 


N4 




2.065(3) A 


2.112(4) A 


Ni - 


N5 




2.062(3) A 


2.070(4) A 


O/S 


- Ni 


- Nl 


88.44(11) 


86.98(15) 


Nl - 


Ni - 


N2 


79.02(13) 


78.04(17) 


N2 - 


Ni - 


N3 


94.62(12) 


98.04(12) 


N3 - 


Ni - 


N4 


79.26(11) 


79.20(15) 


N4 - 


Ni - 


N5 


90.28(12) 


92.46(15) 


N5 - 


Ni - 


O/S 


83.24(11) 


83.30(14) 



perimental methods like SQUID magnetometry, EPR 
and torque magnetometry are possible. Torque magne- 
tometry has been proven to be a very valuable tool in 
this field. HSIl 



II. EXPERIMENTAL 



A. Preparation and Crystal Structures 



The tetranuclear Ni(II) cluster [Ni 4 L|] • 4CH 2 C1 2 
with L b = C5H4N-CON-CN4-C2H5 was prepared as de- 
scribed in Ref. [HJ. The isostructural compound [Ni 4 Lg] 
• 4CH 2 C1 2 with L c = C 5 H 4 N-CSN-CN 4 -C 2 H 5 was syn- 
thesized with a method analogous to that used for the 



[Ni 4 L 



g j complex. 
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FIG. 1: (a) Sketch of both ligands L 6 and L c . (b) Structural representation of [N14LI] and (c) [Ni4Lg] (view along the 
crystallographic SVaxis, H atoms omitted). 



sketched in Fig. 2] They only differ by one position: The 
C=0 group in L h is replaced by a C=S group in L c . 

The crystal structures were determined by X-ray struc- 
ture analysis of single crystals. jToL Both compounds 
crystallize in the space group 74(1) /a. They exhibit 
crystallographic S4 molecular symmetry with the four 
nickel centers forming almost regular squares (Fig. [TJl. 
The molecular S4 symmetry axes are perpendicular to 
the planes of the molecules defined by the Ni centers, 
and necessarily coincide with the magnetic z-axes of the 
molecules. The shape of the crystals is quadratic bipyra- 
midal. The molecular magnetic z-axes are thus parallel 
to the S4 symmetry axis of the crystals. 

The eight ligands in a molecule coordinate the nickel 
centers in two different ways (Fig. A set of four 
ligands links two nickel centers each and builds up the 
[Ni 4 L 4 ] 4+ square-like cores. The second set of four lig- 
ands coordinates the nickel centers at the corners of the 
square cores. Each nickel center is surrounded by six 
donor atoms, five N and one O for [N14LI] and five N 
and one S for [Ni4L§], respectively, forming slightly dis- 
torted octahedral coordination spheres. 

Although the sulfur donors are significantly larger than 
the oxygen donors, the structures of the two complexes 
are remarkably similar. This is evident from a careful in- 
spection of Figs.^b) and^fc). The most notable struc- 
tural difference arises for the CN4-C2H5 groups of the 
corner-ligands. Their orientations differ slightly in the 
two compounds. This is plausible since these groups are 
not coordinated to nickel centers. 

However, the geometry of the nickel coordination 
sphere as well as the structure of the ligand linking two 
nickel ions are almost not affected by the replacement 
of oxygen with sulfur (Table HJ). The Ni-Ni next neigh- 
bor distance is 5.567(5) A for [Ni 4 L|] and 5.560(5) A for 
[NljXg]. Several further distances and bond angles are 
listed in Table [I] for both compounds. 

These structural elements are most relevant for the 
magnetic properties, i.e. ligand-field and superexchange 



interactions. As the two molecules are almost struc- 
turally identical, the potentially different magnetic prop- 
erties of [N14LI] and [Ni4Lg] should be controlled predom- 
inantly by the different electronic properties of oxygen 
and sulfur. 



B. Magnetization Measurements 

For magnetization measurements, a single crystal was 
selected by light microscopy in the mother liquor. To 
avoid decomposition, the crystal was transferred directly 
from the solution into Apiezon grease and mounted on 
a plastic straw. The weight of the crystals was typi- 
cally 100 fig. The background signal of grease and sample 
holder was negligible compared to the signal of the crys- 
tals. Magnetic moment was measured with a MPMS-7 
SQUID magnetometer from Quantum Design. The tem- 
perature range was 1.8- 300 K, the maximum magnetic 
field 5.5 T. The measurements were performed for mag- 
netic fields parallel and perpendicular to the z-axis of the 
crystal. Due to the quadratic bipyramidal shape of the 
single crystals, none of the crystal planes is parallel to 
the magnetic z-axis. Therefore, a proper orientation of 
the sample on the sample holder was difficult. The ori- 
entation accuracy was about 10°. Two samples of each 
compound were investigated. 



C. Magnetic Torque Measurements 

Torque measurements of single crystal samples were 
performed with an appropriately designed silicon can- 
tilever torque sensor which will be described in detail in 
the next chapter. As for magnetization measurements, a 
crystal was selected from the mother liquor, immediately 
covered with grease, and mounted on the torque sensor. 
The typical sample weight was less than 10 fig, but could 
be determined only within an error of 50% due to the 
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FIG. 2: Sketch of the silicon torquemeter described in sec- 
tion lHIl The typical orientation of a quadratic bipyramidially 
shaped crystal sample and its z-axis are shown. For adjusting 
the angle O, the device could be rotated in-situ around the 
magnetic y-axis. 

unknown amount of grease covering the crystals. There- 
fore, the number of molecules in a crystal, X , is known 
only roughly and a calibration of the torque signal was 
not feasible. The torquemeter was mounted on a rotating 
plate which allowed an in-situ orientation of the crystal 
with an accuracy of 0.3°. The sample was mounted on 
the sensor with its z-axis perpendicular to the rotation 
axis, as indicated in Fig. [21 The positioning accuracy was 
better than 3°, which is much better than for the mag- 
netization measurements. Figure [21 defines the angle 
between magnetic field and magnetic z-axis of the crys- 
tal. The torquemeter provides a resolution of 1CP 11 Nm. 
It was inserted into a 15/17 T superconducting cryomag- 
net system. The temperature was adjusted by a variable 
temperature insert. Typically, torque measurements ver- 
sus applied magnetic field were performed at 15 different 
angles in a range of 220°. The lowest temperature was 
1.8 K. For several samples, the low- field range was inves- 
tigated in detail, i.e. at 60 different angles in a range of 
180°. 5 samples of each compound were investigated. 

III. SILICON CANTILEVER TORQUEMETER 

A schematic drawing of the cantilever device is pre- 
sented in Fig. |21 The cantilever is mounted on a sub- 
strate of crystalline quartz glass. It is made from a very 
pure crystalline silicon wafer (resistivity > 3 kVLcm) with 
a (100) surface orientation and thickness of 300 /im. The 
use of crystalline silicon guarantees excellent mechanical 
as well as non-magnetical properties. As shown in Fig. El 
the single crystal silicon cantilever consists of parts at the 
original thickness of 300 /im and strongly thinned parts. 
This has been achieved by different masks on the two 
sides of the wafer. The masks were formed by 1 /xm thick 
Si02 layers grown on both sides of the silicon wafer which 
were patterned by standard photolithography and etch- 
ing with buffered fluorine acid. The silicon cantilever 
itself was etched by hot KOH to a thickness of 10-30 /im. 
A 200 nm gold layer was evaporated on the bottom side 
of the cantilever. Together with gold pads structured on 
the quartz substrate it forms a capacitor and a reference 
capacitor (which is placed at the right hand section of 



the device in Fig. - Finally, the cantilever was glued on 
the quartz substrate. 

The torque of the sample causes a deflection of the can- 
tilever which is detected by a change of the capacitance 
AC. For readout, the two capacitors were connected 
to a ratio transformer forming an ac bridge. With this 
setup a sensitivity of AC /Co = 10~ 7 is readily obtained. 
[l3|,[l4| Co denotes the zero-field capacitance (1 pF). The 
properties of the cantilever torquemeter can be mod- 
eled as follows: The deflection is Ad = (3/2)t/(D c L c ), 
where D c is the spring constant and L c the length of the 
cantilever. The change of the capacitance is given by 
AC/Co ~ Ad/d (l + Ad/d ), where do denotes the dis- 
tance of the capacitor plates. This relation shows nonlin- 
ear behavior. However, for all measurements presented in 
this work nonlinearity was less than 2% and could be ne- 
glected. With Uq being a characteristic of the ac bridge, 
its output voltage can be expressed as U — Uo(AC/Co)- 
Altogether, one obtains U = Kt(1 + aKr) with the cali- 
bration constant K and the nonlinearity a. If required, K 
and a may be obtained from an explicit calibration which 
can be done quite easily in many ways. [Til fl5l Il6l| As 
we were not able to determine the weight of the samples 
accurately, calibration of the device was not necessary. 



IV. THEORY 

The appropriate Hamiltonian for molecular spin clus- 
ters consisting of Ni(II) centers is given by |7L fl7| 

h = - y^jijSi ■ Sj+ 

S, • D« s • S, + fi B J2 • g* • B (1) 

i i 

with Sj = 1 and the following standard terms: a 
Heisenberg term modeling isotropic next-neighbor ex- 
change interactions, a zero-field- splitting (ZFS) term due 
to ligand-field interactions, and the Zeeman term. Due 
to the small coupling constants in [N14LI] and [M4LU, 
anisotropic and biquadratic coupling terms can be ne- 
glected. Dipole-dipole interactions are negligible due to 
the large distances of the spin centers. Cross-coupling 
superexchange terms hardly exist as corresponding cou- 
pling paths are not present. Due to the molecular S4 sym- 
metry of the complexes, magnetic anisotropy is strictly 
uniaxial and a simplified Hamiltonian is obtained: 

H = Si ■ S j + Dj2(Sl ~ 2/3) + 

i<j i 

VB9xy(S x B x + SyB y ) + fi B g z S z B z . (2) 

Thus, four magnetic parameters are sufficient to de- 
scribe the properties of the N14 squares correctly: the 
coupling constant J, the ZFS parameter D and the two 
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FIG. 3: (a) Magnetic susceptibility times temperature vs. 
temperature and (b) magnetic moment vs. magnetic field at 
T = 1.8 K of a [Ni 4 L g ] crystal for two different orientations of 
magnetic field. The solid lines represent best fits based on 
eq. @. 



g-factors g xy and g, 
rameterization g = 



In the following we will use the pa- 

for 



'( 2 9%y + S*)/3 and Ag = g z -y xy 
the g-factors. As we will see, J ~ 0.9 K and D w —2.5 K. 
Thus neither the strong exchange limit (\D/J\ <C 1) nor 
the Ising limit (\D/J\ 3> 1) is valid. Therefore, a full 
matrix diagonalization has to be performed. Due to the 
small dimension of the Hilbert space of 3 4 = 81, this 
can be done on a commercial PC. Calculation time could 
be reduced by a factor of 12 by taking into account a 
Ci v spin permutational symmetry of Hamiltonian eq. • 

m 



V. MAGNETIZATION MEASUREMENTS: 
RESULTS AND ANALYSIS 

Magnetic susceptibility and magnetization curves for 
[Ni 4 Lg] with magnetic field parallel and perpendicular to 
the magnetic z-axis are shown in Fig. [3] Fits based on 
eq. (|2Jl revealed the following parameters: J = 0.40(4) K, 
D = -2.80(5) K, g = 2.33(9) and Ag = 0.01(1) for 
[Ni 4 L|] and J = 0.40(4) K, D = -1.80(8) K, g = 2.33(9) 
and Ag = 0.01(1) for [Ni 4 L|]. It should be noted, that 
the given errors reflect statistical uncertainties only and 
do not include systematical errors. Detailed simulations 
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FIG. 4: Torque measurements vs. magnetic field of a [Ni 4 L|] 
crystal sample at several angles O and T = 1.8 K (circles). 
The curves are shifted for clarity. The solid lines represent fits 
using Hamiltonian eq. The arrows indicate the inflection 
points of the curves. 



showed that the low temperature behavior of the magne- 
tization is very sensitive to a misalignment of the crystal. 
Already a misalignment of 5° , which is well within experi- 
mental uncertainty (see section lTTBTl , leads to notably dif- 
ferent parameters. Therefore, the results from the mag- 
netization measurements will be regarded as guidelines 
and will not be discussed further. Nevertheless, a trend 
is indicated: The coupling constants are roughly identical 
in both compounds, whereas the ZFS parameters differ 
significantly. 



VI. TORQUE MAGNETOMETRY: RESULTS 
AND ANALYSIS 



The torque measurements were analyzed by fitting 
Hamiltonian eq. to the data. As the torque signal was 
uncalibrated, the number X of molecules in a sample had 
to be considered also as a free parameter, in addition to 
the parameters J, D, g, and Ag. In our first attempt to 
determine the magnetic parameters, we fitted the com- 
plete angular dependence t(B, 9) of a sample with X, J, 
D and Ag varying simultaneously, g was fixed to reason- 
able values between 2.1 and 2.3. [7l ll9t 20] A typical data 
set with corresponding fit is shown in Fig.0] Concerning 



the parameter R=V] (r s 



07 E- 



, , which es- 



timates the quality of the least square fits, we obtained 
excellent results. Unfortunately, this approach revealed 
rather large variations of the parameters for different 
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TABLE II: Magnetic parameters for the examined molecules. 





J (K) 


D (K) 


A g 


[NuLg] 


0.9(1) 


-2.7(1) 


0.01(1) 


[Ni 4 m 


0.9(1) 


-2.0(1) 


0.01(1) 



samples. However, fixing one additional parameter be- 
sides g yielded stable fitting results. Thus, either J, D 
or Ag should be determined in a different way. As we will 
show, this can be done without additional experimental 
data. The improved strategy rests on the fact, that a fit 
to the whole data set does not consider that the param- 
eters affect individual parts of the torque curves selec- 
tively. A careful study of numerical simulations based on 
Hamiltonian eq. © provides valuable information con- 
cerning this topic. 

Figure |3{b) shows that the magnetization curves sat- 
urate at fields of about 5T. This is easily understood as 
calculations show that the ground state is well separated 
from the exited states at 5 T (AE>5K). The torque mea- 
surements exhibit similar saturation (Fig.QJ. In this field 
regime, the torque curves are featureless and the magni- 
tude is controlled by X, Ag, and D simultaneously. Here 
the curves are over parameterized. 

For a detailed examination of the low-field part, it is 
useful to analyze the behavior of the inflection points 
B; p (0) of the torque curves. In this way, the parameter 
X is eliminated as the value of Bi p (0) is independent of 
the magnitude of the curves. Figure 01 shows that the 
inflection points vary slightly with O. The calculated 
angular dependence of for different values of J and 
D is presented in Fig. |SJ It turns out that the offset 
of the oscillation is shifted downwards with increasing J 
(and g, not shown here). On the other hand, the offset is 
not influenced by D and, as further simulations showed, 
by Ag. In contrast, the amplitude of the oscillation is 
controlled by D exclusively. So it is possible to extract 
the values of D and J (for given g) from analyzing the 
inflection points of the measured torque curves. As g 
should be between 2.1 and 2.3, J may be obtained with 
an accuracy of 10%. 

For two samples of every compound, we determined 
B, p (0) from measurements at 60 different angles. Fitting 
this data (Fig. HJJ) revealed D and J very precisely. To 
obtain the remaining parameter Ag, the whole torque 
data set t(B, O) was fitted with D and g fixed. As J 
was not fixed in this procedure, it could be determined 
again and compared to the value extracted from Bi p (0). 
This provided prove for the consistency of our analysis. 
The parameters for the two compounds are summarized 
in Table iHl 



VII. DISCUSSION 

The value of J observed for powder samples of [NULg] 
in Ref.^Jis almost a factor of two larger than that deter- 
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FIG. 6: Bi P vs. angle of measured torque curves for (a) 
[NuLg] and (b) [NuLg]. Solid lines represent fits using eq. 
For [NuLg], data sets for two temperatures have been used 
simultaneously in the analysis. 



mined here. We ascribe this to structural changes of the 
cluster resulting from loss of CH2CI2 molecules upon dry- 
ing of crystals. This effect has been observed for other 
molecules, which decompose rapidly under exposure to 
air, too. 0, The different values for J obtained from 
single crystal SQUID measurements as compared to the 
torque results is explained by misalignment errors as dis- 
cussed in section Ivl 

Within the error ranges, the value of Ag is consistent 
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FIG. 7: Sketch of the ligand Lg bridging two nickel ions, the 
free metal orbitals of the two linked Ni centers, and the two 
obvious coupling paths. 



with predictions of ligand-ficld theory: From Ag = 2D/X 
and A = — 250 K appropriate for Ni(II) centers one 
obtains A.g^O.02. 

The above analysis demonstrated that the coupling 
constants are identical within experimental accuracy in 
both compounds, whereas the ZFS parameters differ sig- 
nificantly. The geometry of the Ni coordination spheres 
are essentially identical for the two complexes. Thus, 
different values of D should be ascribed to different elec- 
tronic environments of the Ni centers. In particular, the 
different donor capabilities of oxygen and sulfur should 
clearly affect the ligand-fields and thus the ZFS parame- 
ters. 

The results for J are more puzzling. The special ge- 
ometrical arrangement of the two coordination pockets 
leads to an othogonality of the metal orbitals. To point 
this out, Fig.[7|shows a sketch of the ligand L b linking two 
Ni centers and the relevant metal orbitals (hypothetical 
orbitals of a free Ni atom) . The magnetic orbitals should 
extend along the coupling path between neighboring Ni 
ions as indicated by the gray background in Fig. [7| The 
coupling path actually may be split in two sub-paths, 
path 1 along the Ni-N-C-O-Ni chain and path 2 along 
the Ni-N-C-N-Ni chain. This suggests an explanation for 
the observed ferromagnetic couplings: When path 1 and 
path 2 contribute equally to the magnetic coupling, the 
overlap of the magnetic orbitals will be zero for parity 
reasons. Then, magnetic coupling would be ferromag- 
netic since an antiferromagnetic contribution, which is 
proportional to the overlap, cancels out. [2l| 

However, path 1 of the coupling path is quite differ- 



ent for [N14LI] and [Ni4L§], as demonstrated by the ob- 
served different values for D. Thus, one would expect 
that also the coupling constants are clearly different, in 
contrast to experimental observation. This suggests that 
the magnetic coupling is predominantly controlled by 
path 2. But then, an antiferromagnetic coupling should 
be expected which generally dominates if the overlap is 
non-zero. [211 ] Obviously, the ferromagnetic coupling in 
[NLjLg] and [Ni4L§] is the result of a subtle balance be- 
tween ferromagnetic and antiferromagnetic contributions 
which is not easily reconciled with present knowledge. 



VIII. CONCLUSION 

The two M4 compounds studied in this work are of 
interest for three reasons: (i) Ferromagnetic coupling 
in molecular spin systems is rather rare. Concerning 
polynuclear Ni complexes, only few systems are known 
to date. [2^] (ii) The high crystallographic symmetry re- 
duces the number of magnetic parameters significantly, 
allowing very accurate determination of the magnetic pa- 
rameters, (iii) Two species with minimal geometrical dif- 
ferences of one system facilitates an isolation of possible 
links between ligand-field splitting or magnetic coupling, 
respectively, and electro-structural properties. 

We showed that J and D can be determined very ac- 
curately by torque magnetometry in combination with 
sophisticated analysis also for systems where J and D 
are on the same order of magnitude. This extends re- 
cent applications of torque magnetometry to molecular 
nanomagnets. 000 

The differences of the anisotropy parameter D have 
been ascribed to local differences in the electronic envi- 
ronment of the Ni centers. For the coupling constant J 
some mechanisms have been suggested but no final ex- 
planation could be given. Ab-initio calculations would 
be of great help to provide comprehensive explanations 
for the origin and the strength of the coupling. [23| 
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